i 


AD-A035  706  NAVAL  POSTGRADUATE  SCHCCl.  PONTERfY  CALIF 

AN  examination  of  TURP"LF*  1 r.S«;iPATlON  IN  ThP 
OeC  76  6 « KARCH 

UNCLASSIFIED 


END 

DATE 

FILMED 

3^77 


MICRCX-OMY  HLSOlUIK'N  UbI  < MARI 
NA’l''N-l  hllkMu  <.l  .'AM.AHii  , • 


MCUMTV  CLAMtriCATIOM  Or  THIS  (WMn  Om« 


WFOKT  OOCUMENTATtOH  FACE 


•BAD  MnHUCnONS 
'OKB  COMPI.BTING  FORM 


Naval  Postgraduate  School 
Monterey,  California  93940 


II.  COMmOLLIMS  O^FlCI  NMII  AMO  ADOUCM 

Naval  Postgraduate  School 
Monterey,  California  93940 


14.  MAMiTawMa  AaiNev  mami  a AOONtsvif  «mw«m 

Naval  Postgraduate  School 
Monterey,  California  93940 


TT^af^TI  *"7 


MCURITV  CLAM.  (1  i 

Unclassified 


Approved  for  public  release;  distribution  unlimited. 


IT.  MtTmSUTMM  VTATIMCMT  d B» 


la.  tUPFLIMiMTAIIV  MOT 


Shipboard  measurements  of  temperature  and  velocity  fluctuations  were 
performed  to  determine  optical  propagation  properties  of  the  marine 
boundary  layer.  Erolrl cal  expressions  describing  the  temperature  struc- 
ture parameter ,^C^  In  terms  of  the  Richardson  number,  R1,  overland  were 
used  to  evaluatr-om  obtained  for  open  ocean  conditions.  Profiles  of 
the  dissipation  of  turbulent  kinetic  energy,^  with  respect  to  height, 

Z,  and  with  respect  to  a stability  parameter,  Z/L,  were  examined  for 


DO  I jSTm  M71  taiTIMIDF  IMOVWW 

1)  9/9  I 


I 


lAmi 


tfccuwTY  ccowyic^Tiow  or  tm*  »>a»fw»«w  f>»««  tmi—mt 


OS^vTiA 

In  general,  there  was  little  correlation  between  the  measured/Cr/ 
and  the  stability  parameter.  However,  the  distribution  of  c with'-  ' 
height  for  both  the  stable  and  unstable  cases  showed  little  deviation 
from  that  expected  for  near  neutral  conditions. 


An  Examination  of  Turbulent  Dissipation 
In  the  Marine  Boundary  Layer 


Author 


Approved  by: 


by 


George  William  Karch 
Lieutenant,  United  States  Navy 
B.S.,  United  States  Naval  Acaden\y,  1970 


Submitted  In  partial  fulfillment  of  the 
requirements  for  the  degree  of 


PIASTER  OF  SCIENCE  IN  METEOROLOGY  AND  OCEANOGRAPHY 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1976 


Thesis  AdvUor 


Second  Reader 


^ ~ chairman , bepartment  o 


of  Meteoro'logy 
bean  of  Science  and  Engineering 


3 


ABSTRACT 


Shipboard  measurements  of  temperature  and  velocity  fluctuations  were 

perfonned  to  determine  optical  propagation  properties  of  the  marine 

boundary  layer.  Empirical  expressions  describing  the  temperature  struc- 
2 

ture  parameter,  Cj  , In  terms  of  the  Richardson  number,'  R1,  overland  were 
used  to  evaluate  data  obtained  for  open  ocean  conditions.  Profiles  of 
the  dissipation  of  turbulent  kinetic  energy,  e,  with  respect  to  height, 

Z,  and  with  respect  to  a stability  parameter,  Z/L,  were  examined  for 
open  ocean  conditions. 

In  general,  there  was  little  correlation  between  the  measured 
and  the  stability  parameter.  However,  the  distribution  of  e with  height 
for  both  the  stable  and  unstable  cases  showed  little  deviation  from  that 
expected  for  near  neutral  conditions. 


4 


TABLE  OF  CONTENTS 


I.  INTRODUCTION 11 

II.  THEORETICAL  BACKGROUND - 13 

A.  GENERAL 13 

B.  STABILITY  CONSIDERATIONS  20 

III.  DATA  COLLECTION 23 

A.  PLATFORM  AND  LOCATION  23 

B.  INSTRUMENTATION  AND  PROCESSING  EQUIPMENT  23 

IV.  ANALYSES  PROCEDURES  33 

A.  ANALYSES  OF  MEAN  DATA 33 

B.  ANALYSES  OF  FLUCTUATION  DATA 34 

C.  ANALYSES  OF  VARIANCE  DATA - 43 

D.  HOT-WIRE  CALIBRATION  46 

V.  RESULTS 49 

VI.  CONCLUSIONS 59 

LIST  OF  REFERENCES  70 

INITIAL  DISTRIBUTION  LIST  71 


il 


j 

f 

/ 

'i 

! 


LIST  OF  FIGURES 

1.  Effects  of  unstable  and  stable  conditions  on 

height  profiles  ofO 19 

2.  The  dependence  of  Richardson  number  on  stability  22 

3.  Comparison  of  dimensionless  wind  shear  observa- 
tions with  Interpolation  formulae  22 

4.  Optical  ranges  and  R/V  Acania  position  (normal) 

In  the  Monterey  Bay  .......  24 

5.  Sensor  location  on  board  R/V  Acania  — - — 25 

6.  Thornwalte  Associates  cup  anemometer  wind  pro- 
file register  system  ... 27 

7.  TSI  model  1210  probes  and  wind  vanes  ---  29 

8.  Hewlett  Packard  model  HP-2850  temperature 

sensitive  quartz  crystal  probes  30 

9.  Aspirated  shelter  30 

10.  Dunmore  type  lithium  chloride  humidity  sensor  — - - 30 

11.  Mean  wind  profile  — . — - - 35 

12.  Potential  temperature  profiles  36 

13.  Virtual  potential  temperature  profiles  ---....  37 

14.  Velocity  spectrum  39 

15.  Temperature  spectrum  - --  --  --  --  --  --  --  --  --  - 40 

16.  Spectrum  calibration  plot  - --  --  --  --  --  --  --  --  41 

17.  Sample  strip  chart  section  - 45 

18.  In-si  tu  calibration  plot  - - 43 

19.  Uyngaard  et  a1  (1971)  prediction  curve  - --  --  --  --  --  50 

20.  Observed  prediction  curve  - 51 


7 


21.  Erroneous  profile  >.55 

It  22.  Garrett's  results  -.-56 

23.  Observed  results  -------57 

24.  e versus  logZ ....58 

25.  e versus  logZ------ -.-58 


I 


!k 


i 


LIST  OF  SYMBOLS 


» 


Refractive  Index  structure  function  parameter 
Temperature  index  structure  function  parameter 
Rate  of  dissipation  of  turbulent  kinetic  energy 
Temporal  frequency 
Acceleration  due  to  gravity 
Wave  number 

Monin-Obukhov  stability  length 

Outer  scale,  lower  limit  of  inertial  subrange 

Inner  scale,  upper  limit  of  inertial  subrange 

Atmospheric  pressure 

Separation  distance 

Refraction  index 

Richardson  number 

Spectral  density,  i»f  or  k 

Scaling  temperature 

Mean  horizontal  wind  speed 

Friction  velocity 

Potential  temperature 

Virtual  potential  temperature 

Height 

Stability  parameter 


9 


ACKNOWLEDGEMENTS 


I would  like  to  give  special  recognition  and  thanks  to  Or.  Kenneth 
L.  Davidson  for  his  expert  guidance,  encouragement  and  support  during 
completion  of  this  study.  Many  thanks  also  go  to  Dr.  Thomas  Houlihan 
for  his  technical  assistance  and  advice.  Much  appreciation  is  extended 
to  both  Dr.  Gordon  Schacher  and  to  Dr.  Chris  Fairall  for  their  invaluable 
help  in  the  analysis  of  the  data  collected  for  this  project.  I thank 
Mr.  Steve  Rinard  for  his  careful  analysis  of  the  computer  programs 
developed  for  this  thesis. 


10 


I.  INTRODUCTION 


The  advent  of  more  complex  and  complete  numerical  models  for  analyz- 
ing and  predicting  large  scale  atmospheric  motion  has  also  led  to  better 
specifications  of  the  boundary  layer.  The  boundary  layer  Is  that  region 
In  the  first  kilometer  over  the  sea  which  Is  defined  principally  by 
turbulent  transfer  of  momentum  and  heat. 

Optical  propagation  through  the  atmosphere  Is  affected  by  the  re- 
fractive nature  of  the  medium.  In  addition  to  the  regular  variation  of 
atmospheric  refractive  index  with  height,  there  exist  small  inhomogenei- 
ties In  the  refractive  Index  associated  with  fluctuations  In  the  tempera- 
ture and  velocity  of  the  air.  These  fluctuations  cause  random  phase 
and  amplitude  distortions  In  propagating  wave  fronts  and  thus  degrade 
spatial  and  temporal  coherence  in  the  transmission.  The  magnitude  of 
the  effects  places  limitations  on  optical  system  performance  and  must  be 
Included  In  design  and  execution  considerations. 

Descriptions  of  the  small  scale  fluctuations  which  affect  optical 
propagation  have  not  been  as  complete  nor  In  the  quantity  for  the  over- 
' water  regime  as  for  the  overland  regime. 

Initial  experimental  efforts  to  verify  turbulence  theory  predictions 
were  conducted  over  land.  Measurements  of  wind  speed,  temperature,  and 
humidity  In  those  Investigations  were  made  from  a stable  platform  with 
auxiliary  Instrumentation  well  protected  from  external  Influences. 

Under  these  conditions  relatively  accurate  measurements  of  wind  speed, 
temperature  and  humidity  were  performed  under  a variety  of  weather  condi- 
tions. 


II 


The  marine  environment  on  the  other  hand  presents  a number  of  prob- 
lems, especially  In  accessibility,  platform  stability,  and  sensor  pro- 
tection. Overwater  descriptions  are  necessary,  even  though  considerable 
progress  has  been  made  In  overland  Investigations.  The  necessity  exists 
because  of  the  Increasing  evidence  of  the  Influence  on  atmospheric 
motions  by  features  unique  to  the  oceanic  environment. 

It  Is  the  purpose  of  this  paper  to  compare  overland  and  overwater 
results,  also,  to  evaluate  the  assumption  that  near  neutral  conditions 
apply  for  both  stable  and  unstable  conditions  when  evaluating  the  varia- 
tion of  turbulent  kinetic  energy  dissipation  with  height. 
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II.  THEORETICAL  BACKGROUND 


A.  GENERAL 

Turbulence  properties  of  Interest  In  optical  propagation  are  those 
which  could  be  used  to  describe  the  Intensity  of  fluctuations  In  the 
refractive  Index  and  that  which  could  be  used  to  describe  the  size  or 
scale  of  the  refractive  Index  Inhomogeneities.  The  former  has  been  re- 
lated empirically  to  scintillation,  beam  wander  and  beam  spread,  and 
the  latter  along  with  the  former  to  the  Image  resolution  which  Is  Im- 
portant with  respect  to  seeing  conditions.  Similarity  expressions  re- 
lating mean  conditions  to  properties  of  the  turbulent  regimes  and  hence 
the  temperature  structure  function  and  the  dissipation  rate  of  turbulent 
kinetic  energy  have  evolved  through  a series  of  hypotheses  and  asso- 
ciated measurements. 

An  Important  consideration  within  turbulence  theory  Is  whether  the 
turbulence  Is  anisotropic  or  Isotropic.  In  general,  anisotropy  exists 
within  large  scale  eddies  and  Isotropy  exists  within  small  scale  eddies. 
The  universal  formulae  considered  will  be  valid  only  for  Isotropic 
turbulence. 

Anisotropic  turbulence  Is  very  empirical.  Every  time  the  boundary 
changes,  the  turbulence  changes.  However  at  high  wave  numbers,  which 
corresponds  to  small  scale  eddies,  the  turbulence  should  be  Independent 
of  the  boundaries.  Kolmogorov  (1941)  postulated  that  at  high  wave 
numbers  turbulence  should  be  Isotropic  even  though  It  Is  generated 
from  and  embedded  In  anisotropic  eddies  of  low  wave  numbers.  If 
Reynolds  numbers  are  high  enough,  turbulence  will  adjust  through  Inertial 
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transfer  and  viscous  dissipation  until  a statistical  equilibrium.  Inde- 
pendent of  Initial  conditions.  Is  attained.  This  equilibrium  consists 
of  energy  transfer  Into  the  high  wave  number  range  that  Is  equal  to 
energy  lost  through  viscous  dissipation. 

On  the  basis  of  the  Isotropic  nature  of  small  scale  fluctuations, 
only  one  parameter  Is  necessary  to  describe  the  Intensity  of  the  refrac- 
tive Index  fluctuation  over  many  scales.  It  Is  the  refractive  Index 

2 

structure  function  parameter,  , defined  as 

- [<n(x)  - n(x+r)>]^/r^'^^  (1) 

where  n(x)  and  n(x+r)  are  refractive  Indices  at  two  points  on  a 
line  oriented  normal  to  the  mean  wind  direction  and  separated  by  the 
distance  r which  Is  less  than  the  outer  scale,  L^,  the  lower  end 
of  the  Inertial  subrange,  and  greater  than  the  Inner  scale,  the 
smallest  scale  of  naturally  occurring  turbulence. 

A parallel  expression,  which  defines  the  temperature  structure 
function  parameter  Is 


• [<T'(x)  - T'(x+r)>]^/r^^^  (2) 

where  T'(x)  and  T'(x-«-r}  are  temperature  fluctuations  at  two 
points  separated  by  the  distance  r . The  refractive  Index  Is  deter- 
rnlnted  primarily  by  density  fluctuations  and  can,  therefore,  be 

related  to  temperature  fluctuations,  neglecting  himidity. 

2 2 

is  related  to  Cj  by  the  following 

• C79.X10’®  P/T^]2  (3) 

U 


2 2 

Fortunately,  both  and  Cj  are  readily  measurable  by  optical 

and  meteorological  means,  respectively. ?>n  alternate  relationship  for 
2 

Cj.  , which  involves  measurement  of  the  rates  of  dissipation  of  turbu- 
lent kinetic  energy,  e,  and  temperature  variance,  Xi  is 

(4) 

where  B is  an  empirical  constant  with  a value  of  3.2.  This  last  form 
enables  indirect  estimates  of  to  be  made  from  mean  conditions 
since  c and  x easily  related  to  boundary  layer  fluxes  and, 
eventually,  profiles  if  steady  horizontally  homogeneous  conditions  exist. 

Expressions  which  relate  Cj  to  mean  properties  of  the  boundary 
layer  are  desirable  because  the  small  scale  measurements  are  impractical 
to  obtain  in  most  operational  or  tactical  regimes  and  situations. 

Since  turbulence  is  nearly  synonymous  with  temperature  fluctuations, 
it  is  ultimately  desirable  to  describe  mean  thermal  stratification  in 
terms  of  atmospheric  bulk  stability  parameters  such  as  the  Monin-Obukhov 
length  (L).  In  this  regard,  measurements  of  both  atmospheric  mean 
profiles  and  fluxes  are  required  for  a complete  determination  of  atmos- 
pheric transmission  behavior.  Profile  (ae^/az  and  aO/az)  and  bound- 
ary flux  (U*  and  T*)  parameters  appear  in  the  following  expressions  for 
R1  and  L 


gaSy/az 

e{aO/az)^ 

(5a) 

(5b) 

IS 


/ 


The  folloMing  similarity  predictions  for  the  dependence  of  e and 

X on  momentum  and  heat  fluxes  and  height  were  considered  by  Wyngaard, 

2 

et  a1.  (1971)  In  deriving  an  empirical  expression  for  estimating 
from  mean  stability  parameters 

zl/uj  • fyil/l)  (6a) 

XZ/T*U*«  f2(Z/L)  (6b) 

2 

The  form  of  the  empirical  expression  for  Cy  Is  obtained  by  direct 
substitution  of  equations  (6a  and  b)  Into  equation  (4)  yielding 

C^2  . T*V2/^f3(Z/L)  (7) 

Furthermore  since  (Z/L)  and  R1  can  be  functionally  related, 
Businger,  et  a1.  (1971),  a parallel  dependence  on  R1  can  be  obtained, 
viz.. 


- Z^/^3S/3z)^(|.3(R1)  (8) 

The  functions  f-j  • ^2  , f3  , and  43  In  equations  (6),  (7),  and 
(8)  are  empirical  and  formulated  on  the  basis  of  observations  of  both 
temperature  fluctuations,  momentum  and  heat  fluxes,  and  mean  gradients 
of  wind  speed,  temperature,  and  humidity. 

The  relation  expressed  by  equation  (8)  provides  a desired  depen- 
dence  of  Cj  on  more  readily  measured  mean  stability  (Z,  39/ 3z,  and 
R1).  The  forms  of  f3  and  ^3  and  the  data  obtained  from  the  exten- 
sive AFCRL  study  of  turbulence  structure  over  a flat,  unobstructed 


Kansas  plain  are  presented  In  Uyngaard.  As  will  be  shown  later,  avail- 
able marine  data  does  not  appear  to  agree  for  as  well  as  expected 
with  the  overland  predictions.  I.e.  Hughes  (1976).  The  overland  pre- 
dictions for  f^  and  f2  In  equation  (6),  were  as  follows 

fjCZ/L)  - 4.9(l-7(Z/L))’2/^  Z/L  < 0 

fjCZ/L)  • 4.9(1+2.75(Z/L))  , Z/L  > 0 (9) 

Small  scale  velocity  fluctuation  properties  are  of  Interest  In 
optical  propagation  because  Image  resolution  has  been  empirically  re- 
lated to  the  Innerscale,  which  Is  defined 

Aq  • (10) 

where  y Is  the  kinematic  molecular  viscosity  and  e Is  the  dissipa- 
tion  rate  of  turbulent  kinetic  energy  similar  to  . e can  be  ob- 
tained from  either  one-dimensional  velocity  variance  spectral  estimates 
In  the  Inertial  subrange  or  from  velocity  structure  function  estimates. 

e and  can  be  functionally  related  to  mean  profile  and  flux 
estimates  (U«.  and  R1  or  Z/L)  on  the  basis  of  the  empirical  expres- 
sion presented  previously  (equation  6a).  For  the  purpose  of  examining 
overwater  e results  from  different  stability  conditions  described  by 
Richardson  numbers  (R1),  equation  (6a)  can  be  expressed  as  follows 

e ■ U*^/Kz  ♦,(R1)  (11) 

Since  Z/L  and  R1  are  functionally  related. 
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Overwater  e values  can  also  be  evaluated  on  the  basis  of  U« 
values  estimated  from  them.  Such  estimates  can  be  evaluated  by  compar- 
ing them  with  U«  estimates  computed  from  mean  wind  profile  measure- 
ments. Both  estimates  utilize  the  following  expression  for  the 
mean  wind  gradient 

30/3Z  ■ U^kz  ♦^(R1)  (12) 

where  <|>^(R1)  Is  defined  In  equation  (11).  ((i^(R1)  Is  equal  to  1.0 
under  neutral  conditions.  Is  assumed  to  be  a constant  with  height 
In  this  expression  so  It  Is  restricted  to  the  constant  surface  layer. 
Integration  of  equation  (12)  for  near  neutral  conditions  yields 

0 ■ U*/k  An  Z/Zq  (13) 

I 

i 

I where  Is  assumed  equal  to  zero,  Z^  Is  the  roughness  parameter. 

The  logarithmic  profile  suggested  by  equation  (13)  Is  altered  by 
j stability,  shown  In  Figure  1.  Z^  can  be  eliminated  from  equation  (13) 

by  selecting  mean  winds  at  two  appropriate  levels  so  that 

I 

( 

I U*  - k(02  - 0i)/An(Z2/Z^)  (14) 

I 

It  Is  Important  to  note  that  this  expression  relating  U«  to  mean  wind 
values  (O2  and  0^)  at  two  levels  (Z2  and  Z^)  Is  applicable  only  for 
near  neutral  conditions  (R1>  0). 

In  near  neutral  conditions,  turbulent  kinetic  energy  production  Is 
assumed  to  be  equal  to  the  rate  of  molecular  dissipation  of  turbulent 
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kinetic  energy  and  the  following  relation  Is 
valid 

e » OO/az)  (15) 

combining  equations  (12)  and  (15),  assuming  neutral  conditions  0^(R1) 

« 1.0,  and  solving  for  U*  , yields 

U*  » (ekZ)^/^  (16) 

Thus,  under  neutral  conditions,  the  friction  velocity  (U«)  can  be  esti- 
mated from  either  mean  wind  profiles  using  equation  (14)  or  from 
fluctuation  data  (Involving  turbulent  energy  dissipation)  using  equa- 
tion (16). 

From  determined  values  of  , It  Is  also  possible  to  calculate  a 
momentum  drag  coefficient,  C^,  corresponding  to  a given  height  In  the 
surface  layer 

(17) 

Numerous  studies  have  been  conducted  to. determine  a representative 
value  of  for  10  meters,  Cardone  (196d).  Check  calculations  of 
CiQ  employing  U«  determined  from  both  mean  wind  profiles  and  dissi- 
pation rates  are  performed  In  this  study. 

B.  STABILin  CONSIDERATIONS 

Observational  experiments  by  Businger  et  a1.  (1971)  yielded  a 
definite  relationship  between  the  Richardson  number,  R1,  equation  (5a), 
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and  the  Honin-Obukhov  length,  L,  equation  (5b).  Figure  2 from  Businger, 
et  a1.  (1971)  llllustrates  this  relationship  quite  well.  The  following 
expressions  are  approximations  of  the  relationship  between  Z/L  and  R1, 
proposed  by  Oyer  and  Hicks  (1970)  and  Webb  (1970)  foi  unstable  and 
stable  conditions,  respectively, 

Z/L  - R1 
Z/L  • R1/1-aR1 

where  a Is  an  empirically  derived  constant  equal  to  5.0. 

It  Is  Important  to  note  that  Businger,  et  a1.  (1971)  observed  the 
Richardson  number  to  approach  a critical  value  of  0.21  as  Z/L 
approached  « suggesting  that  as  stability  Increases,  the  flow 
becomes  essentially  non- turbulent.  This  Implies  that  for  the  surface 
layer  the  effect  of  mechanical  turbulence  becomes  negligible. 

The  effect  of  hydrostatic  stability  on  surface  layer  mechanical 
turbulence  Is  further  Illustrated  In  Figure  3 from  Businger,  et  a1. 
(1971).  Recalling  equation  (12)  In  light  of  Figure  3,  It  becomes 
apparent  that  the  function  <^m(Z/L)  valres  relatively  little  with 
respect  to  Z/L  for  unstable  conditions  approaching  a value  of  0.5. 
However,  the  rapid  Increase  of  <^m(Z/L)  as  the  atmosphere  becomes 
hydrostatically  stable  Indicates  that  the  near  neutral  assumptions 
are  Invoked  under  only  slightly  stable  conditions. 


(18) 

(19) 


21 


1 


Figure  2.  The  dependence 
of  Richardson  number  on 
stability. 


Figure  3.  Comparison  of 
dimensionless  wind  shear 
observations  with  Inter* 
polatlon  formulae. 
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III.  DATA  COLLECTION 


A.  PLATFORM  AND  LOCATION 

All  observations  were  made  aboard  the  R/V  Acania  while  anchored  in 
about  30  fathoms  of  water  off  Pt.  Pinos  in  Monterey  Bay.  Open  ocean 
differs  from  land  in  the  effects  of  wave  action  on  turbulence*  in  the 
nature  of  aerosols  and  fog,  and  humidity  fluctuations.  These  conditions 
could  best  be  obtained  far  at  sea;  however  the  cost  of  such  activities 
made  it  desirable  to  work  near  land  and  the  optical  propagation  experi- 
ments could  be  coincident.  Pt.  Pinos  and  other  locations  along  the 
Monterey  Bay  shoreline  provided  very  nearly  the  ideal  situations  for 
both  the  optical  and  meteorological  experiments.  Figure  (4). 

B.  INSTRUMENTATION  AND  PROCESSING  EQUIPMENT 

Measurements  were  made  at  four  levels.  Figure  5 shows  positions 
of  the  towers  which  support  the  instrumentation,  length  of  cable  runs, 
and  position  of  processing  and  recording  equipment. 

The  mean  wind  measurements  were  obtained  with  Thornwaite  Associates 
cup  anemometer  wind  profile  register  systems,  model  number  104.  Three 
cups,  plastic  cones  reinforced  with  aluminum  frames,  are  mounted  on 
stainless  steel  tubes  attached  to  the  main  shaft  at  120®  intervals.  When 
the  cups  rotate  due  to  the  wind  a slotted  shaft  in  the  anemometer  serves 
as  a shutter  between  a light  source  and  a photoelectric  cell.  The  three 
cup  anemometers  have  the  characteristics  of  low  starting  speeds  with  a 
small  amount  of  internal  friction  which  aids  in  checking  inertial  over- 
shoot. 
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Figure  5.  Sensor  location  on  board  R/V  Acania 


Velocity  fluctuation  measurements  were  obtained  using  TSI  model  1210 
probes  with  sliding  support  shields  and  tungsten  wires.  The  shields  per- 
mitted Isolation  of  the  sensing  area  for  the  determination  of  the  undis- 
turbed velocity,  V^,  voltage  reading  before  and  after  each  experiment. 

The  platinum  coated  tungsten  w1  re  was  small  enough  to  resolve  the  viscous 
dissipation  scale  without  making  corrections  for  wire  length.  The  wire 
had  plated  ends  for  Isolating  the  sensing  area  and  thereby  minimizing 
flow  disturbance.  Figure  (6) 

Electronics  associated  with  each  probe  were  a TSI  model  1054B 
linearized  anemometer  and  a TSI  model  1056  variable  decade  module.  The 
anemometer  had  a linear  frequency  response  from  (dead  calm)  0 KHz  to 
10  KHz  and  the  variable  decade  module  operated  with  a 0-60  ohm  range. 

Temperature  fluctuations  were  measured  using  similar  size  sensors. 

A platinum  wire  was  used  In  place  of  a tungsten  wire.  The  temperature 
fluctuation  system  was  designed  for  a resolution  of  0.001°C  at  frequen- 
cies up  to  1 KHz. 

The  temperature  fluctuations  were  measured  using  a bridge  developed 
by  GTE  Sylvania,  the  GTE  Sylvania  model  140.  The  system  was  slightly 
modified  for  use  aboard  the  R/V  Acanla. 

The  baseband  portion  of  this  system  Is  basically  a balanced  wheat- 
stone  bridge  excited  by  a 3 KHz  signal  with  a synchronous  detector  on 
the  output.  Segments  of  a very  small  diameter  platinum  wire  serve  as 
temperature  sensors  In  opposite  arms  of  the  bridge.  The  resistance 
temperature  coefficients  result  In  an  output  from  the  bridge  which  Is 
proportional  to  the  temperature  difference  between  the  two  probes.  The 
sensor  wire  Is  0.5  centimeters  long  and  2.5x10**^  centimeters  In  diameter. 
This  extremely  small  mass  allows  a response  to  temperature  variations 
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Figure  6.  Thomwalte  Associates  cup  anemometer  wind 
profile  register  system. 


of  up  to  1 KHz,  while  electronic  amplification  allows  temperature  differ- 
ences as  small  as  0.004°C  to  be  observed. 

Both  wind  and  temperature  fluctuation  data  were  recorded  on  a San- 
born model  3950  fourteen  channel  tape  recorder.  Real  time  readout  on 
an  eight  channel  chart  recorder,  brush  model  240,  was  used  to  check  the 
quality  of  the  signals  coming  from  the  sensors.  The  chart  data  was  also 
utilized  In  temperature  variance  analyses.. 

Hewlett  Packard  model  HP-2850  temperature  sensitive  quartz  crystal 
probes  were  used  to  obtain  mean  temperature  at  the  sea  surface  and  four 
tower  levels.  Figure  (7)  The  RF  signal  from  the  probes  and  a reference 
oscillator  were  mixed  In  a Hewlett  Packard  HP-2801 A readout  unit  to  pro- 
duce a beat  frequency  whose  signature  can  be  analyzed  to  within  0.001 
°C/Hz.  Each  sensor  received  pre-experiment  calibration  against  a 
platinum  resistance  wire  thermometer  In  a temperature  controlled  circulat- 
ing water  bath.  It  was  checked  over  the  range  of  expected  temperatures. 
The  accuracy  In  achieving  a 0.005“C  correction  factor  was  a constant 
for  each  probe.  The  tower  mounted  sensors  were  housed  In  aspirated 
shelters  to  eliminate  radiation  effects.  Figure  (8) 

Mean  relative  humidity  Information  was  obtained  using  Dunmore  type 
lithium  chloride  sensors.  Figure  (9)  This  sensor  was  also  placed  In 
the  aspirated  shelter.  The  basic  principle  of  operation  of  this  sensor 
Is  resistance  change  In  an  electrolytic  solution  generating  a reference 
voltage  variation  proportional  to  relative  humidity  change.  Automatic 
temperature  compensation  in  the  Instruments  meet  the  following  specifi- 
cations for  relative  humidity,  ± 3S  humidity  below  90S  and  + 4S  humidity 
above  90S.  Sensor  calibration  was  accomplished  by  a comparative  method 
using  a saturated  saline  solution  In  a closed  container. 
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Figure  9.  Aspirated  shelter 


Figure  10.  Ounmore  type 
lithium  chloride  humidity 
sensor. 


Figure  8.  Hewlett  Packard  model  HP-2850  temperature  sensitive 
quartz  crystal  probes. 


Sensor  placement  required  exceptionally  long  cable  runs.  Here  de- 
crease In  system  frequency  response  due  to  their  length  had  little 
effect  In  the  frequency  band  of  Interest.  Adjustments  were  made  In 
the  bridges  for  resistance  and  capacitance  of  the  wire  length. 

Data  logging  during  the  experiments  was  accomplished  using  the  NPS 
developed  MIDAS  (Mlcroprogrammable  Integrated  Data  Acquisition  System). 
This  ml cro-processor  based  data  acquisition  system  utilizes  an  Intel 
8008  central  processor  to  control  the  sampling,  averaging,  and  record- 
ing of  mean  meteorological  data.  All  software  Is  written  In  PL/M  to 
facilitate  the  writing  of  a self-documenting  program. 

The  operator  Is  Interfaced  with  the  system  via  teletype  for  full 
duplex  Input/output  communications  and  program  control  over  the  sample 
start  time  and  the  number  of  samples  to  be  averaged  before  outputting. 
The  operator  may  also  alter  the  present  preset  sample  list  by  adding 
or  deleting  various  sensors  as  they  come  on  line  or  become  Inoperable 
or  Inconsistent.  Once  Initiated,  the  system  Is  fully  automated  to 
sample  the  tailored  list  of  sensors  every  30  seconds  and  periodically 
print  output  values  averaged  over  the  selected  Interval  of  from  one 
minute  to  one  hour.  In  this  study  a 10  minute  averaging  Interval  was 
selected. 

Output  values  were  printed  out  on  the  teletype  In  a columnized  for- 
mat with  the  time  of  print  as  a leader.  The  teletype  has  a paper  punch 
Incorporated  which  may  be  activated  by  the  operator  to  produce  a data 
copy  concurrent  with  the  printout.  This  paper  punch  was  also  utilized 
to  produce  data  cards  for  use  with  the  IBM  360  digital  computer.  A 
magnetic  cassette  tape  recorder  has  been  Integrated  Into  the  system 
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as  a third  data  output  device.  This  cassette  can  be  Interfaced  to  an 
HP  9830  portable  computer  so  that  profile  and  gradient  flux  estimates 
can  be  performed  on  board  automatically  using  the  basic  programming 
capability  of  the  HP  9830. 
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IV.  ANALYSES  PROCEDURES 


4 


A.  ANALYSES  OF  MEAN  DATA 

Mean  wind  data  for  ten  minute  Intervals  were  available  from  the 
Naval  Postgraduate  School  developed  MIDAS  (Mlcroprogramnable  Integrated 
Data  Acquisition  System)  output.  MIDAS  read  values  at  ten  second  In- 
tervals and  averaged  them  for  a ten  minute  period.  Analyses  were  per- 
formed on  data  recorded  on  teletype  printout  and  also  on  the  cassette 
tape. 

The  data  was  edited  for  gross  errors  or  Inconsistencies  due  to 
known  Instrumentation  malfunctions.  The  criteria  for  retaining  or 
discarding  data  from  Individual  levels  or  for  entire  ten  minute  Inter- 
val depended  on  performance  check  results  obtained  during  measurements, 
obvious  Inconsistencies  between  levels  and  sequential  times,  and  the 
availability  of  fluctuation  data.  After  mean  data  were  obtained  for 
ten  minute  Intervals  and  passed  preliminary  editing  for  obvious  erro- 
neous values.  It  was  processed  for  three  applications,  I.e.,  mean  wind 
0,  mean  temperature  T,  and  mean  humidity  q.  These  values  were  plotted 
on  2-cycle  logarithmic  paper.  Since  D,  T,  and  q are  expected  to  vary 
logarithmically  with  height,  a best  fit  straight  line  should  describe 
the  distribution  of  data  points.  Best  fit  lines  were  used  to  estimate 
gradient  parameters  and  the  Richardson  number  (equation  5a,  Figures  10 
and  12).  In  general,  the  procedures  were  subjective  and  In  some  In- 
stances, best  fit  lines  with  different  slopes  could  be  visualized  for 
a given  set  of  data  points.  Therefore,  one  criterion  used  was  not  to 
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give  a single  data  point  too  much  influence  in  determining  the  best  fit 
line.  Consequently*  the  line  drawn  represented  a most  probable  posi- 
tion between  data  points  without  necessarily  coinciding  with  any  point 
(Figure  12). 

8.  AfJALYSES  OF  FLUCTUATION  DATA 

2 

Spectral  analyses  to  obtain  Cj  and  e values  were  performed  on 
fluctuating  velocity  and  temperature  data  obtained  with  single  wires 
oriented  in  the  vertical  normal  to  the  mean  wind.  These  data  represent 
a point  measurement  and  yield  temporal  or  time  descriptions  of  the 
fluctuations.  Twenty-one  minute  segments  of  data  originally  recorded 
on  magnetic  tape  were  recorded  into  the  EMR-Schl umber  model  1510  digital 
spectrum  analyzer.  Procedures  for  converting  spectral  values,  obtained  ' 
with  the  analog  spectral  analyzer,  to  engineering  units  and  for  obtain- 
ing turbulence  parameters  from  the  spectra  are  described  in  the  follow- 
ing paragraphs. 

A necessary  procedure  was  to  scale  the  spectral  plots  in  order  to 
relate  rms  input  voltages  to  power  spectral  densities  (PSD),  variance 
per  unit  frequency.  To  obtain  PSD  levels,  corresponding  to  rms  voltage 
input,  calibrated  scales  charts  had  to  be  constructed. 

For  purposes  of  the  X-Y  plot  format  of  the  analyzer  output,  the  rms 
voltages  were  converted  to  y ■ (voltage)  units  and  a graduated 

scale  was  constructed  enabling  the  logarithm  of  rms  volts  to  be  inter- 
polated from  the  spectral  plots.  The  vertical  scale,  Y,  values  from 
the  plots  were  adjusted  for  each  spectrum  as  a function  of  input  and 
spectral  gains.  The  values  were  then  converted  to  PSD  levels  for  use 
in  calculating  turbulence  parameters.  The  relation  used  was 

(cal.  level  y , 

S(f)  {PSD  units)  . i.5  KnA.{dtr  « <’«  > <“> 
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where  bandwidth  Is  the  frequency  range  divided  by  the  number  of 
channels 

ou  . Freq  range  * 

Number  channel s 

Ok 

» 256  Hz/256  - 1 Hz 

The  calibration  level  volts  rms  equals  the  voltage  at  Y»0.  (1  volt  mms 
for  a 3.16  volt  Input  setting). 

Amplitude  scaling  calibrations  were  accomplished  using  externally 
generated  white  noise.  Signals  with  a 1 volt  rms  over  a frequency 
range  of  0.1  Hz  to  1000  Hz,  giving  a PSD  of  10’^  v^/Hz.  Setting  a 
3.16V  (10  dB)  Input  gain  on  the  spectrum  analyzer  Insures  that  1 volt 
Input  corresponds  to  Y>0.  An  example  of  such  a calibration  plot  Is 
shown  In  Figure  (15). 

2 

The  turbulence  parameters  and  e were  obtained  from  the 
variance  spectra  on  the  basis  of  the  universal  formulae  for  the  Inertial 
subrange  In  wave  number  space.  These  expressions  predict  a -5/3  slope 
for  variance  spectra  of  both  velocity  and  temperature  when  plotted  In 
log-log  fonnat.  Figures  14  and  15  are  typical  spectra  considered  In 
analyses.  The  velocity  variance  spectra  have  consistent  -5/3  slope. 
However,  temperature  spectra  often  exhibited  slopes  slightly  different 
from  the  expected  value  of  -5/3,  This  feature  of  temperature  spectra 
has  been  observed  by  others  and  the  existence  of  "cold  spikes"  In  the 
temperature  traces  has  been  postulated  as  a probable  cause,  Friehe 
(1976). 

Assuming  -5/3  slopes  for  the  variance  specta,  the  Intercept  of  the 
slope  with  the  1 Hz  frequency  line  was  the  power  spectral  density  value 

used  in  computing  the  parameter  of  Interest,  e and  . 
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The  measured  PSD  value  obtained  from  velocity  spectra  was  converted 
to  a spectral  density  In  engineering  units  by  the  equation 

s„(f)  - ■ PSO 

■ (cra/s/volt)^  • volt^/Hz  » (cm/s)^/Hz 

(21) 

where  C Is  the  hot-wire  calibration  factor.  Similarly,  the  equation 
used  to  convert  PSD  values  obtained  from  temperature  spectra  to 
dimensional  values  was 

S^(f)  • • PSD 

■ (®C/volt)^  • volt^/Hz  « •Ahz  (22) 

where  Is  the  calibration  factor  for  the  temperature  system. 

Since  velocity  and  temperature  fluctuations  were  measured  at  a fixed 
point  In  the  flow,  the  resultant  spectral  estimates  correspond  to 
"temporal"  frequencies,  f,  as  Indicated  by  the  notation  S^(f)  and 
Sj(f)  In  the  two  preceding  equations.  To  obtain  e and  values, 
temporal  (f)  and  space  (k)  scales  must  be  related  In  order  to  use 
equations  (21)  and  (22).  This  Is  accomplished  by  Taylor's  (1938) 

"frozen  turbulence"  hypothesis,  that  Is 

k - 2irf/0  (23) 

where  0 Is  the  moan  wind  at  the  moasurement  level.  The  term  "frozen" 
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Implies  the  turbulence  pattern  remains  unchanged  as  It  sweeps  past 
the  sensor  probe. 

The  following  Identles  are  those  used  to  relate  temporal  spectral 
estimates  to  spatial  spectral  estimates  * 


f.Sy(f)  = k*Sy(k)  - c2/3  ,^-2/3 

(24) 

f.s^(f)  = k*S^(k)  » Cg  k'^/^ 

(25) 

where  Sy(f)  and  Sj(f)  are  spectral  density  values  with  units  of 

2 2 

(m/s)  /Hz  and  "C  /Hz  respectively.  The  above  equations  are  rewritten 

2 

In  the  following  forms  to  obtain  e and  Cj  from  values  of  f,  0, 
Sy(f).  and  S^(f) 


• f • S^(f)}2/2 


(26) 

(27) 


where  k * 2rf/0  and  C-j  * .5  and  C2  * .25  , empirically. 

C.  ANALYSES  OF  PAIRED  TEMPERATURE  SENSOR 
2 

Cj  values  were  also  estimated  by  using  paired  sensors  for  variance 
analyses.  This  analysis  was  based  on  the  expression  for  the  tempera- 
ture structure  function,  equation  (2). 

The  analysis  Is  one  In  which  temperature  differences  between  paired 
resistance  wires  separated  by  a distance  r Is  subjected  to  variance 
analyses  as  Indicated  by  the  term  [T'(x)-T'(x<*‘r)]  on  equation  (2). 
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Voltages  corresponding  to  the  teiiepertture  differences,  AT,  «s 
measured  by  the  paired  sensors  were  recorded  on  analog  magnetic  tape. 
During  experimental  periods,  calibration  voltages  were  also  recorded 
on  the  data  tape  for  reconstruction  of  the  signal.  This  calibration 
signal  was  played  back  at  the  time  of  data  analysis  to  provide  a 
reliable  check  of  the  recording  system  performance  In  both  record 
and  playback  modes. 

The  variance  or  mean  square  of  the  temperature  differences  required 

2 

to  compute  was  obtained  by  analyzing  the  analog  voltages  with  a 
mean  square  voltmeter  (Thermo  Systems  Incorporated,  model  1060)  jslng 
a 30  second  time  constant  for  the  averaging.  The  output  of  the  mean 
square  voltmeter  as  well  as  a record  of  the  original  signal  were  recorded 
on  standard  strip  charts.  A strip  chart  section  appears  In  Figure  17. 

The  record  of  the  original  Instantaneous  signal  provides  an  Important 
means  of  spotting  extraneous  noises.  The  rms  voltmeter-strip  chart  re- 
corder system  was  also  calibrated  using  known  Input  voltages.  At  least 
one  Input  calibration  voltage  was  larger  than  the  largest  expected  signal 
on  the  data  record  and  likewise  one  Input  calibration  voltage  was  lower 
than  the  smallest  expected  signal.  This  procedure  minimized  any  non- 
linearities  In  either  the  rms  voltmeter  of  the  strip  chart  recorder. 

Five  minute  averaged  values  of  the  rms  voltages  were  abstracted  from 
the  strip  chart  and  these  values  are  utilized  In  a program  for  the  IBM 

360  computer.  This  program  converts  the  mean  square  voltages  values  to 

2 

temperature  difference  variances  values,  from  which  Cj  was  computed. 

A final  check  on  the  validity  of  the  measured  data  was  made  by  perform- 
ing spectral  analyses  on  representative  data  samples  from  each  experi- 
mental run.  The  spectrum  analyzer,  described  In  section  B,  was  used  to 


obtain  the  spectra  which  were  subsequently  plotted  and  kept  as  a 
permanent  record.  Any  harmonic  noise  appeared  as  a spike  In  an  other- 
wise reasonably  smooth  curve.  It  should  be  emphasized  that  continuous 
calibration  and  data  checking  procedure  were  essential  steps  through- 
out the  analyses. 

0.  HOT  WIRE  CALIBRATION 

In-si tu  calibrations  of  the  velocity  sensors  were  accomplished 
every  10  minutes  during  the  experiment.  In  this  procedure  recordings 
were  made  of  both  the  cup  anemometer  wind  speed  and  the  corresponding 
hot  wire  voltage  output.  The  sensor  wind  speed  Is  given  by 

■ a + b (28) 

where  V Is  the  hot-wire  voltage  output  and  U Is  the  mean  wind  speed 
for  any  given  level.  The  constants,  a and  b,  are  the  calibration  curve 
slope  and  Intercept  respectively.  The  constants,  a and  b,  are  derived 
from  the  In-si tu  calibration  curve. 

To  convert  the  voltage  PSD  levels  to  velocity  units  requires  a 
calibration  factor  given  by 

e 

U'  ■ Cv'  (29) 

where  C Is  the  calibration  factor  In  m/s  /volt,  v'  Is  the  voltage 
fluctuation,  and  U*  the  velocity  fluctuation.  Differentiating  equa- 
tion (28)  yields 

U*  ■ (4vtf^^/a)  v'  (30) 
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which  when  substituted  Into  equation  (29)  yields 

C » 4vtf ^^/a 

An  exanple  In-situ  calibration  curve  Is  shown  In  Figure  18.  The 
slope  of  the  curve  yields  the  value  for  a In  equation  (31). 


(cm/sec) 


VOLTS* 


Figur*  18.  In-situ  calibration  plot. 
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V.  RESULTS 


A comparison  of  the  observed  overwater  results  and  those  obtained 
by  Wyngaard,  et  al.  (1971)  appear  In  Figures  19,  20,  respectively. 

In  each  case  dots  represent  an  Individual  data  point.  For  the  overwater 
results  mean  values  over  R1  Intervals  of  0.25  were  determined,  plotted 
In  the  center  of  the  Interval,  and  appear  as  small  dots  within  circles. 
Error  bars  delineate  one  standard  deviation  from  the  mean  within  each 
Interval.  The  number  at  the  top  of  each  error  bar  Is  the  number  of 
data  points  used  In  defining  the  mean  In  the  Interval. 

For  the  stable  and  unstable  regime  with  0.5  > R1  >-0.5  there  Is 
some  correlation  with  the  results  of  Wyngaard,  et  al.  However,  for 
R1  <-0.5  there  Is  little  agreement  with  the  difference  In  results  on 
the  order  of  one  order  of  magnitude.  There  Is  some  agreement  with  the 
trend  In  the  sense  that  It  doer  not  change  significantly  with 
Instability. 

Many  of  the  Richardson  numbers  were  not  within  acceptable  ranges, 
for  positive,  and  negative  R1  values.  Results  are  not  Included  for 
positive  R1  values  above  1.0  and  negative  values  lower  than  -3.0.  A 
predicted  critical  positive  R1  value  can  be  derived  from  equation  (19), 
which  gives  a value  of  0.21. 

« 

Scatter  In  the  observed  results  could  be  attributed  to  the  scatter 

2 

In  both  the  measured  Cj  and  3e/3z  values.  Deviation  of  tempera- 
ture spectra  from  the  -5/3  slope  caused  uncertainty  In  the  Cj^  esti- 
mates. The  3d/3z  values  are  dependent  on  profile  estimates  whose  de- 
termination was  at  times  subjective  as  discussed  In  conjunction  with 
Figure  21. 
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Friction  velocity  values,  U*  , were  calculated  in  two  ways.  In 
the  first  method  (flethod  A)  profiles  of  log  (z)  versus  mean  wind,  U , 
were  plotted  to  obtain  gradient  estimate.  An  example  profile  is 
shown. in  Figure  10.  After  the  profiles  were  drawn,  5 and  10  meter 
level  wind  speeds  were  estimated  from  the  profiles  for  each  10  minute 
period  and  these  values  were  used  to  solve  equation  (14)  for  . 

These  values  appear  in  column  A of  Table  I. 

Friction  velocity  values  were  also  calculated  from  fluctuation 
data  (Method  B).  First,  equation  (24)  was  solved  to  obtain  a value 
for  the  turbulent  kinetic  energy  dissipation  rate,  c , at  each  level. 
Then  the  neutral  prediction,  equation  (16),  was  used  to  obtain  a 
at  each  level.  Ideally  U«  at  all  levels  would  be  equal  for  any 
given  period  since  is  assumed  to  be  constant  in  the  surface  layer. 
The  values,  with  about  10%  exceptions,  were  within  acceptable  limits 
of  each  other.  They  were  then  averaged  to  obtain  representative 
values  for  . Results  from  these  computations  appear  in  Table  I 
column  B. 

values  were  also  calculated  using  wind  speed  dependent  C^q 
values  suggested  by  Cardone  (1969)  and  a value  of  U taken  from  pro- 
files, at  a height  of  10  meters,  using 

U,  . (32) 

Cardone' s C.|q  representation  was  for  the  marine  boundary  layer. 

These  values  are  shown  in  Table  I column  C. 

The  results  are  similar  to  those  reached  by  Atkinson  (1976),  i.e., 
showing  some  disagreement  between  U«  values  obtained  by  the  two 
methods,  with  the  fluctuation  data  method  more  closely  agreeing  with 
Cardone* s results. 
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_3 

The  normalized  dissipation  rate  KeZU*  In  relation  to  Z/L  was 

tested  In  the  same  manner  as  Garratt  (1972).  The  results  of  Garratt 
(1972)  and  those  obtained  In  this  study  are  shown  In  Figures  22  and  23, 
respectively.  The  curve  represents  the  function  (^jjj(Z/L)-Z/L  . Similar 
scatter  was  also  observed  In  an  Investigation  by  Busch  and  Panofsky 
(1968). 

An  examination  was  also  made  of  the  variation  of  the  observed  dissi* 
patlon  rate  of  turbulent  kinetic  energy,  e , with  height  for  the  differ- 
ent stability  conditions.  Results  of  this  appear  In  Figures  24  and  25 
for  stable  and  unstable  conditions,  respectively.  The  results  from 
both  stability  regimes  are  compared  to  the  -1  properties  for  neutral 
stability.  The  results  show  that  for  all  stability  conditions  a slope 
I corresponding  to  near  neutral  conditions  was  observed. 

I 


I 
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VI.  CONCLUSIONS 


2 

CorrtUtlon  of  spectrally  derived  Cj  results  with  the  stability 
paraaeter»  R1.  was  poor  for  the  data  obtained  In  this  study  of  the 
Marine  envlroneent.  Verification  evaluations  for  friction  velocity 
results  froM  Mean  (profile)  and  turbulent  (dissipation)  data  were  satis- 
factory. Correlation  of  noriMlIzed  e with  the  stability  parameter 
(Z/L)  showed  the  same  scatter  as  results  of  others.  The  variation  of 
e with  height  exhibited  results  which  support  the  near  neutral  pre- 
dictions for  most  cases. 

It  Is  concluded  that  the  turbulent  field  In  the  marine  boundary 
layer  Is  subjected  to  several  anomalous  effects  which  can  cause  the 
disruption  of  the  Inertial  subrange.  Thus,  empirically,  expressions 
for  describing  overland  transmission  characteristics  will  undoubtedly 
have  to  be  altered  to  provide  operational  determinations  for  utiliza- 
tion In  the  marine  environment. 


I 

t 
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TABLE  II.  Cj^  and  R1  Results. 


Period 

1974.  3-27-1841 

3-27-1902 
3-27-1928 
3-27-2117 
3-28-1902 
3-28-2059 
8-13-1700 
8-13-1721 


Level 

dS/dz 

(M) 

®C/an 

4.19 

0.049 

6.60 

0.031 

7.64 

0.027 

13.89 

0.015 

4.19 

0.047 

6.60 

0.030 

7.64 

0.026 

13.89 

0.014 

4.19 

0.047 

6.60 

0.030 

7.64 

0.026 

13.89 

0.014 

4.19 

0.019 

6.60 

0.012 

7.64 

0.012 

13.89 

0.006 

4.19 

0.046 

6.60 

0.029 

7.64 

0.025 

13.89 

0.014 

4.19 

0.054 

6.60 

0.034 

7.64 

0.030 

13.89 

0.016 

4.19 

-0.036 

6.60 

-0.023 

7.64 

-0.020 

13.89 

-0.011 

4.19 

-0.036 

6.60 

-0.023 

7.64 

-0.020 

13.89 

-0.011 

^10-3 

R1 

1.59 

0.362 

0.94 

0.500 

1.57 

0.500 

- 

1.000 

1.97 

0.249 

0.61 

0.392 

2.02 

0.454 

- 

0.750 

5.40 

0.230 

1.37 

0.362 

2.89 

0.419 

1.20 

0.750 

3.45 

0.073 

6.89 

0.115 

- 

0.113 

1.42 

0.242 

16.56 

1.000 

- 

1.500 

12.07 

1.500 

1.46 

1.500 

56.73 

0.388 

- 

0.500 

14.21 

0.500 

19.25 

1.000 

4.33 

-0.538 

0.42 

-0.847 

0.53 

-0.981 

0.23 

-1.783 

3.76 

-0.665 

0.42 

-1.047 

0.67 

-1.213 

0.27 

-2.204 
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Period 

1974  8-13-1742 

9-18-1643 
9-19-1735 
9-19-1846 
9-19-1900 
9-19-1929 
9-19-2011 
9-19-2115 


Level 

dS/dz 

(M) 

®C/cm 

4.19 

-0.025 

6.60 

-0.016 

7.64 

-0.013 

13.89 

-0.007 

4.19 

0.031 

6.60 

0.020 

7.64 

0.017 

13.89 

0.009 

4.19 

-0.005 

6.60 

-0.003 

7.64 

-0.003 

13.89 

-0.001 

4.19 

-0.005 

6.60 

-0.003 

7.64 

-0.003 

13.89 

-0.001 

4.19 

-0.010 

6.60 

-0.006 

7.64 

-0.005 

13.89 

-0.003 

4.19 

-0.005 

6.60 

-0.003 

7.64 

-0.003 

13.89 

-0.001 

4.19 

-0.010 

6.60 

-0.006 

7.64 

-0.005 

13.89 

-0.003 

4.19 

-0.172 

6.60 

-0.109 

7.64 

-0.094 

13.89 

-0.052 

i^io-3 

l\l 

4.25 

-0.451 

0.45 

-0.711 

0.63 

-0.823 

0.21 

-1.496 

0.29 

0.308 

0.08 

0.486 

- 

0.500 

0.11 

1.000 

0.47 

-0.082 

1.18 

-0.128 

- 

-0.149 

0.56 

-0.270 

5.07 

-0.100 

10.27 

-0.157 

- 

-0.181 

1.10 

-0.330 

0.59 

-0.607 

0.87 

-0.957 

- 

-1.107 

0.91 

-2.013 

0.73 

-0.099 

1.17 

-0.156 

- 

-0.180 

1.14 

-0.328 

0.50 

0.051 

0.81 

0.080 

- 

0.093 

0.55 

0.169 

0.35 

-3.000 

0.69 

-3.000 

- 

-3.000 

0.77 

-3.000 
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Period 

Level 

dS/dz 

c 2 

Ri 

(M) 

®C/cm 

x'lO  ^ 

9-19-2158 

4.19 

-0.010 

0.46 

-0.202 

6.60 

-0.006 

1.12 

-0.319 

7.64 

-0.005 

- 

-0.369 

13.89 

-0.003 

1.10 

-0.671 

9-19-2243 

4.19 

-0.178 

0.32 

-3.000 

6.60 

-0.113 

0.78 

-3.000 

7.64 

-0.098 

- 

-3.000 

13.89 

-0.054 

1.08 

-3.000 

9-20-0014 

4.19 

-0.011 

0.66 

-0.461 

6.60 

-0.007 

2.12 

-0.727 

7.64 

-0.006 

- 

-0.841 

13.89 

-0.003 

1.0 

-1.530 

9-20-0056 

4.19 

-0.013 

0.56 

-0.408 

6.60 

-0.008 

- 

-0.643 

7.64 

-0.007 

* 

-0.744 

13.89 

-0.004 

1.88 

-1.353 

9-20-0140 

4.19 

-0.013 

0.77 

-0.546 

6.60 

-0.008 

1.43 

-0.860 

7.64 

-0.007 

- 

-0.995 

13.89 

-0.004 

0.89 

-1.809 

9-20-0530 

4.19 

-0.003 

0.76 

-0.085 

6.60 

-0.002 

1.77 

-0.134 

7.64 

-0.002 

- 

-0.155 

13.89 

-0.001 

1.68 

-0.282 

9-20-0552 

4.19 

-0.006 

0.91 

-0.952 

6.60 

-0.004 

- 

-1.500 

7.64 

-0.004 

- 

-1.736 

13.89 

-0.002 

0.99 

-3.000 

11-21-1800 

4.19 

-1.626 

3.74 

-3.000 

6.60 

-1.032 

8.85 

-3.000 

7.64 

-0.892 

2.02 

-3.000 

13.89 

-0.490 

5.93 

-3.000 
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Period 

Level 

(M) 

1974 

11-21-1821 

4.19 

6.60 

7.64 

13.89 

11-21-1842 

4.19 

6.60 

7.64 

13.89 

11-21-1903 

4.19 

6.60 

7.64 

13.89 

11-21-2027 

4.19 

6.60 

7.64 

13.89 

11-21-2048 

4.19 

6.60 

7.64 

13.89 

11-21-2109 

4.19 

6.60 

7.64 

13.89 

11-21-2130 

4.19 

6.60 

7.64 

13.89 

1975 

3-27-1503 

4.19 

6.60 

7.64 

13.89 

dS/dz 

*C/cm 

x'lO  ^ 

R1 

0.008 

3.39 

0.027 

0.005 

1.97 

0.043 

0.004 

1.97 

0.050 

0.002 

2.24 

0.091 

0.006 

3.49 

0.033 

0.004 

2.17 

0.051 

0.004 

2.11 

0.060 

0.002 

9.17 

0.108 

0.003 

5.65 

0.009 

0.002 

3.26 

0.014 

0.002 

3.00 

0.017 

0.001 

1.24 

0.030 

0.002 

7.75 

-0.001 

0.001 

1.42 

-0.001 

0.001 

3.61 

-0.001 

0.000 

0.89 

-0.002 

0.002 

8.39 

0.006 

0.001 

9.75 

0.010 

0.001 

3.55 

0.012 

0.000 

1.36 

0.021 

-0.000 

9.71 

-0.025 

-0.000 

7.22 

-0.039 

-0.000 

4.11 

-0.046 

-0.000 

1.19 

-0.083 

0.0 

11.04 

-0.037 

0.0 

7.97 

-0.059 

0.0 

5.90 

-0.068 

0.0 

1.72 

-0.123 

0.006 

10.57 

0.004 

0.004 

6.14 

0.007 

0.004 

3.43 

0.008 

0.002 

- 

0.015 
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Period 


R1 


1975 


Level 

(M) 

3-27-1519  4.19 

6.60 

7.64 

13.89 

3-27-1539  4.19 

6.60 

7.64 

13.89 

3-27-1600  4.19 

6.60 

7.64 

13.89 

3-27-1644  4.19 

6.60 

7.64 

13.89 

3-27-1656  4.19 

6.60 

7.64 

13.89 

3-27-2020  4.19 

6.60 

7.64 

13.89 

3-27-2045  4.19 

6.60 

7.64 

13.89 

3-27-2106  4.19 

6.60 

7.64 

13.89 


d§/dz 

“C/ctn 

. 

x'lO 

0.006 

9.66 

0.004 

5.90 

0.004 

1.65 

0.002 

- 

-0.002 

10.04 

-0.001 

5.86 

-0.001 

2.58 

-0.000 

- 

-0.005 

10.71 

-0.003 

3.42 

-0.003 

2.65 

-0.001 

m 

-0.002 

11.25 

-0.001 

7.46 

-0.001 

3.01 

-0.000 

- 

0.0 

11.31 

0.0 

5.74 

0.0 

2.08 

0.0 

0.041 

0.026 

- 

0.022 

2.76 

0.012 

- 

0.044 

5.17 

0.028 

- 

0.024 

0.20 

0.013 

- 

0.015 

5.45 

0.009 

- 

0.008 

2.01 

0.004 

- 
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« 0.004 
0.006 
0.007 
0.012 

-0.059 

-0.092 

-0.107 

-0.195 

-0.032 

-0.051 

-0.059 

-0.107 

-0.014 

-0.023 

-0.026 

-0.048 

-0.013 

-0.020 

-0.024 

-0.043 

0.016 

0.025 

0.029 

0.053 

0.025 

0.039 

0.045 

0.082 

-0.046 

-0.072 

-0.083 

-0.151 


TABLE  III 
2 

April  Cruise  Cj  and  R1  Results 
Date  Level  1 


Time 

36/ 3Z 

rc/m) 

Cj-^xlO^ 

R1 

4-27-1519 

.135 

4.439 

-0.68 

1529 

.150 

3.288 

-1.15 

1544 

.160 

2.877 

-0.88 

1554 

.160 

2.466 

-0.73 

1625 

.160 

2.634 

-1.42 

1635 

.175 

2.404 

-2.04 

4-28-1433 

.050 

5.178 

-0.75 

1443 

.075 

4.439 

-0.93 

1458 

.050 

7.398 

-0.60 

1517 

.085 

4.809 

-0.80 

1531 

.085 

5.918 

-1.88 

1540 

.050 

5.918 

-20.01 

1613 

.075 

5.548 

-0.94 

1618 

.075 

5.548 

-0.90 

1623 

.105 

4.439 

-1.84 

1633 

.065 

5.918 

-1.98 

1703 

.035 

2.877 

-0.37 

1720 

.125 

3.288 

5.95 

1730 

.925 

3.699 

0.70 

1740 

.020 

3.699 

-0.47 

1909 

.040 

2.774 

1.07 

1919 

.040 

2.959 

4.45 

1929 

.030 

3.144 

5.01 

4-29-1415 

.025 

1.099 

3.34 

1425 

.020 

3.625 

0.50 

2035 

.040 

3.625 

-0.99 

2050 

.040 

3.625 

-1.48 

4-30-1315 

.040 

4.285 

-0.71 

1325 

.040 

5.933 

-1.80 

1335 

.050 

5.273 

-0.61 
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Data  Level  2 


Time 

36/ 3Z 
(“C/m) 

Cy^xlO^ 

Ri 

4-27-1519 

.095 

4.279 

-1.07 

1529 

.105 

2.445 

-1.81 

1544 

.10 

2.649 

-1.38 

1554 

.105 

2.852 

-1.15 

1625 

.105 

1.630 

-2.24 

1635 

.11 

1.936 

-3.22 

4-28-1433 

-1.18 

1443 

.055 

3.260 

-1.46 

1458 

- 

- 

-0.95 

1517 

.060 

4.482 

-1.26 

1531 

.060 

3.871 

-2.96 

1540 

.035 

3.871 

-31.52 

1613 

.050 

3.667 

-1.48 

1618 

.055 

4.279 

-1.42 

1623 

.065 

3.464 

-2.89 

1633 

.015 

4.279 

-3.12 

1703 

.030 

2.649 

-0.58 

1720 

.955 

3.056 

9.37 

1730 

.001 

3.260 

1.11 

1740 

.025 

3.056 

-0.74 

1909 

.01 

2.751 

1.68 

1919 

.015 

2.037 

7.01 

1929 

.005 

- 

7.89 

4-29-1415 

.03 

4.265 

5.26 

1425 

.005 

4.436 

0.79 

2035 

.035 

1.706 

-1.57 

2050 

- 

- 

-2.34 

4-30-1315 

.020 

2.730 

-1.12 

1325 

.020 

2.218 

-2.83 

1335 

.050 

3.071 

-0.96 
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Data  Level  3 


Time 

36/3z 

(*C/m) 

C^^xlO^ 

R1 

4-28-1340 

.045 

3.129 

-1.69 

1350 

.045 

2.869 

-1.37 

1400 

.040 

2.608 

-1.24 

1410 

.040 

2.869 

-0.91 

1423 

.035 

2.347 

-1.21 

1803 

.035 

0.652 

-0.68 

1814 

.035 

0.652 

-0.68 

1827 

.015 

0.782 

-0.61 

1836 

.01 

0.782 

0.27 

1847 

.01 

3.912 

0.34 

1858 

.01 

9.779 

0.34 

4-29-1700 

.065 

2.539 

-7.12 

1711 

.08 

2.537 

-4.7 

1722 

.075 

3.173 

-12.6 

1732 

.08 

3.173 

-5.98 

1744 

.065 

2.539 

-5.98 

4-30-1125 

.05 

3.491 

-2.92 

1137 

.015 

2.856 

0.43 

1232 

.01 

1.904 

0.0 

1242 

.015 

1.523 

0.0 

1252 

.01 

1.142 

0.65 

1302 

.015 

1.269 

0.0 

Data  Level  4 


Time 

30/3Z 

("C/m) 

C^^xlO^ 

Ri 

4-27-1519 

.075 

2.658 

-2.24 

1529 

.075 

1.772 

-3.82 

1544 

.080 

1.772 

-2.90 

1554 

.085 

1.772 

-2.41 

1625 

.080 

1.329 

-4.71 

1635 

.085 

1.476 

-6.77 

4-28-1433 

.015 

4.958 

-2.49 

1443 

.035 

2.953 

-3.07 

1458 

.020 

6.325 

-1.99 

1517 

.040 

4.616 

-2.65 

1531 

.040 

2.067 

-6.22 

1540 

.025 

2.658 

-66.33 

1613 

.030 

2.658 

-3.11 

1618 

.035 

2.953 

-2.99 

1623 

.055 

2.067 

-6.09 

1633 

.045 

. 

-6.57 

1703 

.015 

1.772 

-1.23 

1720 

.075 

2.067 

19.71 

1730 

.020 

2.067 

2.33 

1740 

.010 

1.772 

-1.55 

1909 

.040 

2.137 

3.54 

1919 

.025 

1.476 

14.76 

1929 

.025 

1.476 

16.61 

4-29-1415 

.005 

2.531 

11.06 

1425 

.015 

2.847 

1.65 

2035 

.005 

- 

-3.29 

2050 

.005 

- 

-4.92 

4-30-1315 

.015 

-2.35 

1325 

.015 

- 

-5.96 

1335 

.005 

- 

-2.03 
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TABLE  IV 

April  Cruise  e Results  (e  x 10^) 


Time 

Level  1 

Level  2 

Level  4 

4-28-1330 

3.61 

4.92 

1.62 

1339 

3.97 

2.90 

1.66 

1358 

3.58 

3.67 

1.30 

1757 

4.37 

4.55 

1.93 

1815 

5.83 

4.66 

1.56 

1833 

4.68 

4.95 

1.88 

1851 

4.33 

3.94 

1.87 

4-29-1652 

12.4 

4.00 

6.68 

1710 

8.41 

2.33 

5.83 

1728 

10.4 

3.04 

4.17 

1916 

6.92 

11.6 

2.44 

1934 

5.04 

1.30 

1.80 

2048 

- 

2.21 

4.35 

2203 

17.8 

3.12 

1.15 

2221 

19.4  ' 

3.43 

1.78 

4-30-1013 

4.87 

3.43 

0.84 

1031 

4.82 

0.97 

0.88 

1232 

7.40 

2.09 

1.93 

1250 

3.67 

1.37 

1.37 

i 

t 

1 

I 
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